Relaxation due to random collisions with a many-qudit environment 
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We analyze the dynamics ol a system qudit of dimension /i sequentially interacting with the v- 
dimensional qudits of a chain playing the role of an environment. Each pairwise collision has been 
modeled as a random unitary transformation. The relaxation to equilibrium of the purity of the 
system qudit, averaged over random collisions, is analytically computed by means of a Markov chain 
approach. In particular, we show that the steady state is the one corresponding to the steady state 
for random collisions with a single environment qudit of effective dimension = vfi. Finally, we 
numerically investigate aspects of the entanglement dynamics for qubits {jj, = v — 2) and show that 
random unitary collisions can create multipartite entanglement between the system qudit and the 
qudits of the chain. 

PACS numbers: 03.65.Yz, 03.67.Mn, 03.67.-a 
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I. INTRODUCTION 

The repeated collision model is a simple yet instruc- 
tive model of irreversible quantum dynamics. First in- 
troduced in [l| to analyze the process of thermalization 
and more generally of homogenization 0, 0] in the 
limit of an environment with a large number of degrees of 
freedom, it has been further studied to elucidate various 
aspects of the irreversible dynamics of quantum systems 
in thepresence of environments with few degrees of free- 
dom [a, Q . In most of the literature cited above the inter- 
action between system and environment is due to pair- 
wise elastic collisions modeled by partial swap operators. 
In the present paper we will instead model such collisions 
by random unitary operators 0, H, [13, HHj ■ Such 
choice is due to several reasons. On the one hand, al- 
though an exact modehzation of the system-environment 
interaction is difficult, a good description of the relax- 
ation process can be obtained by a suitable average over 
random interactions. Some examples of such approach 
can be found in Il3 | , where the irreversible dynamics 
of a single and a pair of qubits is analyzed in terms of a 
random interaction with an environment which is itself 
modeled as a random matrix and in 0| , where the average 
dynamics of a single qubit interacting with a very small 
reservoir is described again by random collisions. On 
the other hand random and pseudo-random states and 
their efficient generation by suitable se quen ces of random 
gates [il, [il 0, [il in H, HE m, El,]!! have received 
a considerable attention due to possible applic ations in 
quantum information processing psi . [2^ |27|. |28|| . 

In the present paper we analyze the approach to equi- 
librium of a system qudit interacting with a very large 
ensemble of qudits. The interaction is modeled by a 
sequence of two-qudit random collisions. The paper is 



structured as follows: in the following section we review 
the random collision model and we specialize it to the 
case of collisions described by random unitary operators. 
We then characterize the approach to equilibrium by ana- 
lytically calculating the purity of the system steady state 
and the rate of approach to such state for the specific case 
of colliding qubits. Such analytical analysis is then gen- 
eralized to a system qudit of dimension fi colliding with 
an ensemble of qudits of size v. We then proceed with 
a numerical analysis of the entanglement dynamics for 
colliding qubits 



II. THE RANDOM COLLISION MODEL AND 
THE RELAXATION TO EQUILIBRIUM 

In the random collisions model of irreversible dynamics 
a system qudit interacts with an environment consisting 
of N qudits. Such interaction is modeled by pairwise col- 
lisions between the system qudit an a single environment 
qudit. Each collision is described by a random unitary 
operator. The environment, i.e. N, is assumed to be so 
large that the system never collides twice with the same 
environment qudit. In pictorial terms one can think of 
a single qudit colliding in sequence with the individual 
qudits of a long chain. The overall state of the system 
and environment, after t collisions, is 



(t) 



Uo 



(1) 



where U^j is a random unitary operator acting on the pair 
of qudits 0, j; labels the system qudit and j = \, - ■ ■ , N 
labels the environment qudits. Let us assume that the 
system and environment are in an initial tensor product 
state Q^g'l, = ''r?i?72 ■ ■ - rj^. Since the collision operators 
are random unitaries, the specific states r]i are irrelevant 
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and we can assume, without loss of generality that all the 
environment qudits are in the same initial state rj. 

In order to characterize the relaxation process we first 
consider the decay of the system purity after t collisions. 
We remind the reader that, given a density operator g, 
its purity is defined as 7^ = Tr [gi^] . The purity is a 
decreasing function of the degree of statistical mixture 
of g and, for a qudit of dimension ^ takes values in the 
range - < V < 1, where V = 1 corresponds to pure 
states and V — to the completely unpolarized mixed 
state. Since we are focussing our attention to the purity 
of the system, our model is equivalent to a system qudit 
colliding with a single environment qudit whose state is 
refreshed to its initial state rj after each collision. We will 
show that, after averaging over random unitary collisions, 
1^(1) can be analytically calculated. 



A. Colliding qubits 

Let us first consider the case in which both system and 
environment consist of qubits. The density operator of 
the system and environment qubits can be written as 



QSE 



E ' 



(2) 



where ctq" (ct^^) denotes a Pauli matrix acting on the 
system (environment) qubit, with e {0,x,?/,z} and 
= I. The purity of the overall system and environ- 
ment after t collisions then reads 



and the system's purity is given by 



(3) 



(4) 



Note that Vse is invariant under unitary evolution, i.e., 
the overall system-environment purity is the same before 
and after each collision (of course before the state of the 
environment qubit is refreshed). The constraints 



Tr [gsE] = 1, Tr [g^j^] = Vse 



(5) 



where 



and 



M = 



(2 2 
1^00 7 %x 



(\ 

1^ 



(8) 
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(9) 



Vo 1^ ••• hi 



is a Markov 16 x 16 matrix. The structure of M is clear: it 
leaves Cqq unchanged while all the other components of (? 
are uniformly mixed. The equilibrium state of such chain, 
if the state of the environment qubit is not refreshed after 
each collision, must lie in the subspace spanned by the 
eigenvectors of M with eigenvalue 1. A normalized basis 
corresponding to the unit eigenvalue of M is given by 



(1,0,.. .,0), 



1 



(0,1,...,1). (10) 



The equilibrium state then reads 

V = XqVq -I- X\V\, 



(11) 



where 

xx = {vi,v) 



^00 



1 

16' 



1 



Therefore, 
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-"^ (Qo,ai)#(0,O) 



-iVsE - 1 

16' 



1 ^PsE-1 



16 



15 



240 



4:VoE - 1 
240 



(12) 



Finally, we obtain 



'P^»T.<o = l + ^i'i'PsE~l). (13) 



For pure system-environment states {Vse = 1) we re- 
cover Lubkin's results [291, 



Vl 



(14) 



lead do 



1 



(ao,ai)#(O,0) 



E 



4VsE - 1 
16 



(6) 



It has been shown f2ll| that when two qubit collide 
with a sequence of random U (4) unitaries the ensemble 
averaged coefficients c^(t) evolve according to a Markov 
chain as 



(7) 



It is important to note that such (ensemble averaged) 
equilibrium state is reached after a single random colli- 
sion. In other words, as it should be, the state of the 
system and of the first environment qubit after a single 
collision is a two-qubit random state. 

We now consider the case in which the state of the 
environment qubit is refreshed after each collision. After 
t collision Eq. (fT3|) becomes 



m = l + ^[i'PsE{t-l)-l], (t=l,2,...). (15) 
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Note however that Vse changes when the environment 
qubit is reset to its initial state -q. Just before the t-th 
cohision we have 



VsE{t~l)^V{t-l)V^, 



(16) 



where 7^,, is the purity of the unperturbed environment 
qubit state 77. 

We now solve the equation 

nt) = \ + ^[^V{t-l)rr,-ll (i=l,2,,...). (17) 



The equilibrium value 

V{oo) 



5-2V„ 



is obtained as the solution to the equation 
If we define ^{t) —V{t) — V{oo)^ we obtain 



m = \v,,^{t - 1) = ] m- 



(18) 



(19) 



(20) 



Note that V{po) as well as the decay rate to equilibrium 
depend on P^: the lowest the purity of the environment 
qubits the fastest is the approach to equilibrium and the 
lowest the equilibrium value of the system purity. Al- 
though this has some analogies with the relaxation to 
thermal equilibrium it is worth stressing that the two 
processes have some important differences. In partic- 
ular note that, contrary to the case of thermalization, 
7^(00) ^ Vrf- Furthermore also the dependence on Vri of 
the decay rate is not the one expected for the relaxation 
to equilibrium. This is not surprising as collisions mod- 
eled by random collisions are not elastic. Indeed it has 
been shown in [l[ that the only operator guaranteeing 
thermalization (and in general homogenization ^ ) is the 
partial swap. 

When the chain qubits arc initially in a pure state, we 
have 



V{oo) 



and equation ([^0]) becomes 



(21) 



V(t) 



V{0) 



(22) 

where A — ln(5/2) = 0.916... is the rate of approach to 
equilibrium. Starting from a pure state we obtain V{1) = 
4/5, P(2) = 54/75 = 0.72, ViS) = 258/375 « 0.688, 
in agreement with our numerical data, as shown in Fig.[T] 
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FIG. 1: (color online) Ensemble averaged purity functional 
behavior for /j. = 2,v — 2. The limit value is 'P(cxd) = 2/3. 
The markers are the ensemble averaged values. The curve 
shows the analytic result of Eq. (|22|l 



B. Colliding qudits 

The above Markov chain technique can be generalized 
to analyze the more general case in which a system qu- 
dit of dimension /i collides with a large number of en- 
vironment qudits of dimension v. In this case system- 
environment interactions are modeled by random uni- 
taries drawn from the [/(i)-invariant Haar measure in 
the overall Hilbert space of size L = /ii^. The steps of the 
analysis previously carried on for qubits can be straight- 
forwardly followed to obtain the equilibrium value and 
the decay rate of the purity of colliding qudits. We use 
again representation (f2|), where the cr's are now gener- 
alized Pauli matrices poL [3lj | . The purity of the overall 
system then reads 



ao—O ai —0 

and the system's purity is given by 

OiQ 

The constraints 



(23) 



(24) 



TT[gsE] = l, TT[glE]=VsE (25) 

lead do 

. ,x 1 2 fJ-i^'PsE - 1 



Cool 



(qo,qi)^(0,0) 

Also in this case the vector 

„2 _ /2 



C — (CoO: ' ' ' J C„y) 



(26) 



(27) 



evolves in time according to a Markov chain like in 
Eq. ([7]). We extend Znidaric's conjecture [2l| about the 
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form of the Markov matrix by assuming again that M 
leaves Cqq unchanged while it mixes uniformly all other 
components. M must therefore be a {^jlv) x [piv) matrix 
of the form 



1 



M 





n 1 . . . ____ 



(28) 



The two-dimenisonal eigenspace corresponding to the 
unit eigenvalue of M is spanned by the vector basis 



vq = (1,0, ...,0), vi = 



1 



1 



(0,1,...,!). (29) 



Following the same steps that lead to Eq. (jl2|) we obtain 
that the equilibrium state is 



1 



HvVsE - 1 



^ivVsE ~ 1 



(H2'(/,:.)2[(^^)2_l]'-'(^^)2[(^^)2_l] 



The system purity of such state is 



(30) 



[liuVsE-l). (31) 



If the initial system-environment state is pure {Vse = 1) 
we recover Lubkin's result [29|: 



V = 



^1 + V 



(32) 



This means again that after a single collision the state 
of the system and of the first colliding qudit is a random 
state in the Hilbert space of dimension ^®i> 



V 



12/33 




t 
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FIG. 2: (color online) Ensemble averaged purity functional 
behavior for /i = 4, = 2. The limit value is Vioo) = 12/33. 
The markers are the ensemble averaged values. The curve 
shows the analytic result of Eq. (|37)l 





FIG. 3: From top to bottom, reading from left to right: 
P(1),7'(2),P(3),P(4),P(5),P(6) statistical distributions. In 
the lower side of each figure are shown the statistical average 
(indicated as e) and the standard deviation (indicated as S) 
values. 



If the environment is refreshed after each collision, we 
have 

7'(i + 1) = - + -J-^^^\p-vV{t) - 1]. {t = 0, 1, ...). 



[I ii{{[ivY - 1] 

This equation has steady state 

/i + {[iv) 



Vioo) 



fx{fiiy) + 1 ■ 

If we define ^(i) =V{t) — V{oo), we obtain 

m = <w = «*e(o), 

where 



a = 



Therefore, 



r{t) = e-^* {V{Q) - V{oo)) + V{oo), 



(33) 
(34) 

(35) 

(36) 

(37) 



with A = - In a. In Fig. [21 7^(00) for = 4, = 2 is 
plotted. 
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C. Purity statistics 

Eq. ([M)) suggests that the average system purity is the 
same one would obtain from a partition of a random state 
of a Hilbert space of dimension i.e. Eq. ([34]) coin- 

cides with Lubkin's result if we assume that the system 
qubit has interacted with an environment of effective size 
i^e — [IV. To support this conjecture we have plotted in 
Fig [3] the numerically generated histograms of the statis- 
tical distributions of Vit) for the first 6 collisions for the 
case of colliding qubits. Indeed, as shown in Fig. 21 af- 
ter just six collisions the histogram practically coincides 
with the purity distribution of a system qubit colliding 
with a single qudit of dimension 4. 



e =0.666, 5=0.100 

r r 



FIG. 4: Histograms comparison. Grey shadow: system purity 
distribution after the sixth collision (after which the equilib- 
rium distribution is in practice already reached). Black: pu- 
rity distribution of a system qubit colliding with a single qu- 
dit of dimension = 4, with average = 0.666 and variance 
&r = 0.1000 

An intuitive explanation of such behavior can be given 
in simple terms. After the first collision between system 
and environment qudits the system becomes mixed - even 
if originally it was in a pure state. Following a standard 
mathematical procedure however its state can be purified 
by introducing a fictitious qudit of dimension /i entangled 
with the system qudit. Therefore the whole process can 
be seen as a sequence of collisions among three qudits - 
the third one being the fictitious purification qudit - in 
an overall pure state. There is evidence ^] that in this 
case the overall system will evolve into a pure random 
sate in a Hilbert space of dimension /i^z/. This explains 
the purity statistics of the system qudit. 

We have further numerical evidence to support the 
above analysis: in [sj] Scott and Caves reported an an- 
alytical expression of the average variance of the system 
purity distribution for a random interaction between a 
system of /i degrees of freedom and an environment of v 
degrees of freedom. Such variance turns out to be 



(Aii/ + 3)(Ai// + 2)(/ij/+ 1)2 



(38) 



In our case instead the variance of the system purity dis- 
tribution turns out to be well described by the above 
formula if we make the substitution v — v[i i.e. 



(7?2) - l(Pf = 



1 



(39) 



In the following Table|T]we report the standard deviations 
for different values of and v. In the second column we 
report the value of the Scott & Caves modified standard 
deviations of Eq. ( I39p while in the third column we re- 
port the corresponding numerically computed standard 
deviation values. Indeed a very good agreement between 
the two set of values is clearly seen. 







V 




Scott & Caves 


Collision model 




= 2, 


V 


= 2 


0.1005 


0.1010 




= 2, 


V 


= 3 


0.0769 


0.0767 


M 


= 3, 


u 


= 2 


0.0608 


0.0628 


M 


= 4, 


u 


= 2 


0.0382 


0.0388 




= 2, 


V 


= 4 


0.0618 


0.0632 




= 3, 


V 


= 3 


0.0433 


0.0438 



TABLE I: Comparison between the Scott & Caves modified 
standard deviations (/i, v ^ fj.,Vs = vfi) and numerically cal- 
culated standard deviations 



III. ENTANGLEMENT DYNAMICS 

To further characterize the approach to equilibrium, 
in this section we illustrate some aspects of the entan- 
glement dynamics for qubits, i.e. ^ = v — 2. Since the 
overall state of the system and chain remains pure, the 
entanglement dyanamics is conveniently characterized in 
terms of the so-called tangles [s^ [s^. We remind the 
reader that, given the density operator pij of a bipartite 

is defined as 



system of two qubits, the tangle r^^j 

Ti\j{p) = [max{0, ai — a2 — o-z — 04}]^, 



(40) 



where {0;^} (/c = 1, ... 4) are the square roots of the eigen- 
values (in non-increasing order) of the non-Hermitian op- 
erator = Pij{(Jy ® 0-y)p*^{ay ® (J y) , (J y IS thC tj-V SMIi 

operator and p*j is the complex conjugate of py, in the 
eigenbasis of the ® Uz operator. The concurrence C 



is defined simply as = 



The tangle Ti\j, or 



equivalently the concurrence Cij can be used to quantify 
the entanglement between the pair of qubits i,j for an 
arbitrary reduced density operator pij. In our case the 
overall state of the system and the chain is pure. There- 
fore, the amount of entanglement between qubit i and 
all the remaining qubits can be quantified by the tangle 
Tj|j.gg^ = 4detpi. After t collisions, the tangle ''"0 1 chain 
between the system qubit and the chain conveys the same 
information as the purity 'P{t). Indeed, it is easy to show 
that 



^olchain(^) = 2 - 27^(0. 



(41) 
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The purely multipartite entanglement tm established be- 
tween the system qubit and the qubits of the chain can 
be quantified as 

t 

TMit) = 7-Q|phainW - Xl^olj- (42) 
i=o 

i.e. as the amount of entanglement which cannot be as- 
cribed to purely bipartite entanglement between the sys- 
tem qubit and each individual environment qubit. Due to 
the complexity of the analytical expressions of the tan- 
gles, we resort to numerical simulations. Note that, in 
order to evaluate the tangle '''oichain ^® have to numer- 
ically find the pairwise tangle between the system qubit 
and each qubit of the chain. This implies that we have 
to retain the overall system-chain density matrix; i.e. in 
contrast to the computation of the system purity V , we 
cannot trace over the chain after each collision. In other 
words, as far as the entanglement dynamics is concerned, 
the qubit-chain model is not equivalent to a model in 
which the system qubit collides with a single environ- 
ment qubit whose state is refreshed after each collision. 

The pairwise tangles tqij are shown in Fig.O as a func- 
tion of the number of collisions. A non-zero tangle ro|t' 
between the system qubit ant the i'-th qubit of the chain 
is generated immediately after the t'-th collision and then 
quickly decays for t > t' . The dependence of TQ\t[i) on 
the number t of collisions if fitted satisfactorily by the 
exponential curve To|((t) = 0.166 -I- 0.512 exp(—0.921t), 
suggesting the asymptotic value to|oo(oo) — 0.166 for the 
paiwise tangle generated between the system qubit and 
a qubit of the chain immediately after their collision. 



system qubit and the qubits of the chain as 

TAf((^)« [2-27'M]-To|oo(oo)« i. (43) 

This expectation is confirmed by the numerical data 
shown in Fig. [51 the convergence of TM{t) to its asymp- 
totic value Tm (oo) is well fitted by the exponential decay 
tm(<) - tm{oo) = -1.491 exp(-0.849i), with the fitting 
parameter tm{og) — 0.472 w i. 



0.47 

0.36 

T M 
0.20 




■ 1 2 3 4 5 6 7 8 
t 

FIG. 6: (color online) Multipartite entanglement tm for the 
system qubit as a function of the number t of collisions. The 
solid line represents the exponential fit TKiit) — tm (oo) = 
— 1.491 exp(—0.849t), with the fitting parameter tm{c>o) = 
0.472. 



0.26 

0.20 
0.166 




FIG. 5: (color online) Pairwise tangles To\j as a function of 
the number t of collisions (dashed curves with circles). The 
top circles refer to To\t{t) immediately after the t-th colli- 
sion and are fitted by the solid curve Tomtit) — to|oo(oo) ~ 
0.512 exp(—0.921t), with the fitting parameter ro|oo(oo) ~ 
0.166. 

Our numerical data shown in Fig. [5] also suggest that 
the weight of the terms TQ\j{t) with j ^ t can be ne- 
glected with respect to TQ|j(i). Under this approximation 
and using Eqs. (|^^ and (|^T|l we can estimate the asymp- 
totic multipartite entanglement established between the 



IV. CONCLUSIONS 

We have investigates the dynamics of a system qudit of 
dimension /i sequentially interacting with the qudits, of 
dimension i/, of a chain. Each pairwise collision is mod- 
eled as a random transformation drawn from the Haar 
measure on U{^v). The relaxation to equilibrium, in 
terms of ensemble average over random collisions, in an- 
alytically investigated by means of a Markov chain ap- 
proach. We have shown that the steady state is the one 
corresponding to the steady state for random collisions 
with a single environment qudit of effective dimension 
Ve = i^^-- Furthermore, in contrast to the case of the ho- 
mogeneization process induced by purely elastic partial 
swap collisions random unitary collisions can gener- 
ate multipartite entanglement. 
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